1 0 2 8 8 1 09 68 NX VOVN 


i 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3950 



CHARTS FOR THE ANALYSIS OF FLOW IN A WHIRLING DUCT 

By Robert A. Makofski 

Langley Aeronautical Laboratory 
Langley Field, Va. 


i 



Washington 
May 1957 


t 


AUVL'C 

ULI;„ ,1 (" >r, 


AhVL 


Al 


I 2811 


TECH LIBRARY KAFB, NM 



NATIONAL ADVISORY COMMITTEE FOR AERONAL 


TECH LIBRARY KAru, i« 



TECHNICAL ROTE 3950 

CHARTS FOR THE ANALYSIS OF FLOW IN A WHIRLING DUCT 
Hy Robert A. Makofski 


SUMMARY 


Charts are developed for the analysis of flow in a whirling constant- 
diameter duct. These charts permit the determination of the duct-exit 
Mach number , or the Mach number at any point along the duct, for practical 
ranges of duct-inlet Mach number, duct-tip Mach number, friction factor, 
ratio of ambient temperature to duct-flow total temperature, ratio of 
duct -tip radius to duct hydraulic diameter, and heat conductivity through 
the duct wall and surrounding material. The method of using the charts 
is illustrated through the computation of a sample problem. 


INTRODUCTION 


In certain types of rotor tip -mounted propulsive systems being used 
in helicopters, for example, the pres sure -jet system, the air supply is 
ducted through the rotor blade. In order to obtain a complete analysis 
of such a system, a knowledge of the flow properties in the duct iB 
required. The usual method of analysis requires a time-consuming numer- 
ical integration of a nonlinear differential equation similar to those 
derived in references 1 and 2. The purpose of this paper is to perform 
this integration and to present the results in chart form for a range of 
parameters which include the duct-inlet Mach number, the duct-tip Mach 
number, the duct-inlet temperature, the heat-transfer rate, the friction 
factor, and the ratio of the duct hydraulic diameter to the duct-tip 
radius. The charts are presented for a constant -diameter duct; however, 
possible applications of these charts for ducts of varying diameter are 
discussed. 

The analysis presented herein is limited to steady-state, one- 
dimensional, compres Bible, vIbcous flow and applies only to straight 
ducting. 


SYMBOLS 


A duct cross-sectional area, sq ft 


a 


speed of sound, fps 
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specific heat at constant pressure, ft-lb/slug-°R 

bk 

duct hydraulic diameter, , ft 

Perimeter 

local heat-transfer coefficient based on difference between wall 
temperature and adiabatic wall temperature, ft-lb/ft^-sec- °R 

friction factor, rj i pu^ 

thermal conductivity, ft-lb/ft-sec-°H 

duct Mach number 

duct -tip Mach number, flr e /a a 

heat added to fluid element per unit mass per unit surface 
area, ft-lb/ft^-slug 

heat transferred from fluid element per second, — — 

sec 

thermal resistance between duct wall and outer wall of blade 
per running foot of blade, sec-°R/ft-lb 

radial distance from center of rotation to fluid element, ft 

radial distance from center of rotation to exit of duct being 
considered, ft 

surface area through which heat is transferred per running 
foot, sq ft 

static temperature , °R 

total temperature , °R 

duct-flow velocity, fps 

ratio of fluid -element radius to duct-tip radius, r/r e 
ratio of specific heats 

ratio of duct heat-transfer coefficient to overall heat- 
transfer coefficient (defined in eq. (5b)) 

mass- density, slugs/cu ft 
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t friction shearing stress, lb/sq ft 

Mach number function defined in equation (8b) 
ft rotor angular velocity, radians/sec 

Subscripts : 
a ambient air 

aw adiabatic wall 

e duct-exit station 

i duct-inlet station 

o center line of rotation station 

ow outer wall, that is, wall exposed to ambient air 

w duct wall 


AMItfSIS 


In order to develop an expression for the steady, one-dimensional, 
compressible, viscous flow in a whirling duct, it is necessary to analyze 
all the forces and energies acting on a fluid element. An analysis of 
this type using the momentum and continuity equations yields a relation, 
as given in reference 1, between the duct Mach number and the independent 
flow variables. This relation may be written in the form 



which gives the effect of duct area, friction, whirling, and temperature 
on the duct Mach number. A solution to equation (l) is -possible for 
known values of the independent variables. However, the total -temperature 
variation along the duct is dependent upon Mip and K^Tg/D, and it is, 
therefore, necessary to analyze the energies involved in the flow system. 
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The Energy Equation 

For a steady -flow process, the energy equation is a balance of all 
the energies involved in a system. Therefore, for a fluid element in a 
duct subjected to a centrifugal force in the flow direction, the energy 
equation for constant specific heat is 


dT-p r\ lys 

- Qpr - Q = 
dr D 


( 2 ) 


The heat transferred from the duct may be expressed as 


q ^ JlAp _ rp \ 
H D D pul^ W 


The adiabatic wall temperature T aw is assumed to be equal to the total 
temperature; that is, the recovery factor is assumed to be equal to one. 
The error involved in such an assumption is very small at low duct Mach 
numbers . 


Substituting the preceding expression for the third term of equa- 
tion (2), nondimensionalizing, and using the relation a a 2 = (7 - l)cpT a 

reduces equation (2) to 


*T t /T t 

dx 


, To 


(7 - l)— Mij^x + 
T t 


4h r e fow 

Cppu D \T t 


(3) 


Equation ( 3 ) may be further simplified by use of Reynolds analogy which 
states that 


h = Kf 

Cppu 2 


This relation is fairly accurate at subsonic Mach numbers (refs. 3 and 4). 
Equation ( 3 ) may then be expressed as 


dT t/?t,i 


(7 - 1) 


?a 

^t,i 


Mrp^X + 2 


~*t X&efcv 

T t ,i D 


(4a) 


dx 
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For adiabatic flow in the duct Tv “ T-fc , the total temperature is 
given by 

-in . 1 + ufuz . XJ 2) ( to) 

If the flow is not adiabatic, the ratio of wall temperature to total 
temperature must be determined as a function of x. 

Ratio of wall temperature to total temperature .- An analysis of the 
heat transfer from the duct to the ambient air results in a relation 
between the temperatures of the duct flow, the duct wail, and the ambient 
air. This relation may be expressed as 



Tt 



(5a) 


where k is the ratio of the duct heat-transfer coefficient (gas to duct 
wall) to the overall heat-transfer coefficient (from the gas to the ambient 
air) and is given by 

1 hSrj- 

4=1 + + hSwR (5b) 

b a Sow 

Equation (5a) assumes the aerodynamic heating of the external skin of the 
rotor to be negligible. 

The ratio of the duct heat-transfer coefficient to the overall heat- 
transfer coefficient is determined by the heat-transfer coefficient of the 
gas in the duct, the heat-transfer coefficient between the ambient air and 
external skin of the rotor, and the internal construction of the rotor 
blade. It varies from zero for adiabatic flow to about 0.5 for thin- 
walled ducts in the leading edge of a rotor blade. 

The main difficulty in determining a value of k for a duct in a 
helicopter rotor blade would be in estimating the thermal resistance 
between the i nn er surface of the duct and the outer surface of the rotor 
blade. However, the effect of heat transfer on the change in Mach number 
in the duct is small compared with the effect of whirling and friction. 

It is possible, therefore, to use an approximate value of the thermal 
resistance to compute k and the effects of heat transfer on the Mach 
number variation in the duct. 
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Total -temperature variation In duct .- The variation of total tem- 
perature in the duct may be obtained by substitution of equation (5a) 
into equation (4a) which gives 


(y . iws. ^ . i fii 

T t,i 1 D \ft,i %i/ 


(6) 


For a constant-diameter duct, the solution to equation (6) is of 
the form 


T t 


. 2E 5£i( x - Xi ) (r - ^ “r 2 


T. 


= e 


t,i 




-2k^£(x- X i ) 

; D - 1 + 


o- K f r e 
2k - ir - 


*f*e, 

-2k — - — 


x i e 


— ' 


) 

. o 

- 

T t,i 


1 - e 


-2r^£(x-x ± ) 


(7a) 


where k is a mean value of k for the duct. The term Kj.r e |D is con- 
stant for a constant-diameter duct except for the small effect of the 
temperature variation along the duct on Kj. 

Equation (7a) may be simplified for low rates of heat transfer by 
expanding the exponential terms and keeping only the first two terms of 
the expansion. The resulting equation is 


3_ = 1 + 2LlJ: Jk_ ^2 L.Z _ 2\ + 2 - l\ (x - x^ (7b) 

*t,± ^^ T t,i v 1 D lpb,i y v y 

KKf r e f t . \ 

For values of = 0.007, the value of — 1] computed by equa- 

D / 
tion (7b) is within 4 percent of the value computed by equation (7a)* 

The absolute value of well within 1 percent of the value 

computed by equation (7a)* 
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With, the knowledge of the total- temperature variation along the duct 
as given by equations (6) and (7a.) , it is now possible to write a general 
equation describing the flow in a whirling duct. 


The General Flow Equation 

The general equation describing the flow in a whirling constant- 
diameter duct may be obtained by substituting equation (6) into equa- 
tion (l) and simplifying. Thus, 


dM 2 

dx 


M 2 (l + 


21^1 M 2 ' r 


1 - M c 


4tM 2 


Kf* e 


- (r + 1)% 2 


T„ 


L t,i 


*t,i *b 


x - 


2(1 + 



(8a) 


Equation (8a) is subject to the same assumptions as equations (l) and (6) 
namely, constant specific heat and molecular weight, the validity of the 
Reynolds analogy, and the equality of the adiabatic wall temperature and 
duct total temperature. 

The Mach number at any point in the duct or at the duct exit is 
dependent upon five parameters: M^, Itp, Kfr e ^D, 8X1(1 

Accordingly, any attempt to present in chart form the results of inte- 
grating equation (8a) would involve a large number of calculations and 
plots . 


A change of variable, however, permits results to be obtained in 
terms of only four parameters with satisfactory accuracy. Equation (8a) 
is first rewritten as follows : 


(l - M 2 )dM 2 


d0 


Kf r e/D m 2 ^! + — 1 m 2 ^ K f r e A 1 


4yM 2 - (7 + 1)— — £ 

K f r eA* T t ,i T t 


x - 


2(1 + tM 2 )k (l - 

V T t,i 


dx 


(8b) 
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where 


0 = lo^ 


M 2 


1 + 


a 


_2 


(r+i)/(r-i) 


is the new variable, and 


Mr^ 


is a new parameter. A plot 


Kf r e/-D Kfr e jl) 

of 0 as a function of Mach number is given in figure 1 for y = l.hO. 

The left side of equation (8b) is an exact dif ferential so that its 
integral is not dependent upon the precise variation of M 2 along the 
duct. On the right side, the second (centrifugal) term is normally of 
the order of 5 to 10 times as much as the first or third terms (friction 
and heat transfer, respectively) . Furthermore, the first and third terms 
have opposite signs and, therefore, tend to cancel. Accordingly, the 
variation of M 2 along the duct (which affects only the first and third 

^ — which, in this 


terms on the right side) does not greatly affect 


case, is nearly proportional to 


Mr 2 


Kf.r e /D 


The preceding paragraph dealt with large values of 


Mrp c 


It is 


Kfrg/D* 

of interest to consider the opposite extreme - that of zero rotational 
speed. In this case the second term disappears and only the first and 
third terms remain. At reasonable values of Kfr e jD, there is little 
variation of M 2 along the duct for this case of zero rotational speed 

d0 


so that the integral of 


between the duct inlet and the position 


Kfr e /D 

being considered is relatively small. 

The essential conclusion from the preceding discussion is that it 
should be possible to choose some typical values of Kj>r e /D, to integrate 

d - d 

equation (8b) numerically to get 


Mr]! 2 


%r e /D 


as a function of the four 


parameters , , and K, and to obtain results that are 

within practical accuracy for any other value of Kf r e jD within reason. 
Of course, for large values of Kfr e /D and high duct Mach numbers, this 
reasoning is Invalid. 
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Calculation of Whirling Effect at Low 
Duct Mach Numbers 

As may be seen from figure 1, 0 varies rapidly with Mach number 

at low duct Mach numbers - A small error in the calculation of the duct- 
inlet and duct-exit Mach numbers can cause a significant error in the 
difference between 0 O and 0 e . Therefore, the charts presented in 

this paper are not extended below a duct Mach number of 0.15- For duct 
Mach numbers below this value, it is suggested that the effects of fric- 
tion and whirling be computed separately and added. The Mach number 
variation may then be written 

Me 2 - M± 2 = (^ 2 - Mj 2 )^^ + (Me 2 - Mj 2 )^ (9) 

The first term on the right side of equation ( 9 ) represents the change 
in Mach number of a nonwhirling duct. This value may be computed from 
equation (8b) with M T = 0 or by the method given in reference 5* The 

second term on the right side represents the change in Mach number due 
to whirling and may be computed by integrating the rotation term of 
equation (8b). That is. 


(01 - 0a) 


^ » (r + D % 2 ^ 


T a f 1 T t,i 

^t,i^ X£ T t 


x dx 


( 10 ) 


Since the effect of k is very small at the low Mach numbers, 
equation (*t-b) may be used to express T tjij'Ff Equation (10) then 

becomes 


(01 - 0a) 


= 7 + 1 
Mp 7 - 1 


i°g 6 


1 + 


7 - 1 


T„ 


L t,i 


m t 2 ( 1 - x i 2 ) 


( 11 ) 


The decrease in Mach number due to whirling may be computed by use of 
equation (ll) and figure 1. 

Calculations of the duct-exit Mach number by this procedure show 
good agreement with data obtained by integration of equation (8b) and 
with experimental data. 
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CHARTS FOR ANALYZING THE FLOW 
IW A WHIRLING DUCT 


The charts for analyzing the flow in a whirling constant -diameter 
duct are presented in figures 2 and 5* The curves were computed by 
numerically integrating equation (8b) from Xj_ = 0 to x = 1.0 for 

Kf>r e /D = 0 . 06/. Check points for values of K^r e /D other than 0.06 were 
also calculated. The values of were determined for the various 


values of k. 


T„ 


T. 


t, o 


-, Mq, and 


Kfr e /D 

Hp 2 


W ] 


assumed. 


0 - 

In figure 2, the positive values of — —2 an a function of 

K^r e /D * 


are plotted for a range of values of 


‘•a 


0 to 0-5/ and 


Miji^ 


T. 


t,o 


from 0-95 to 0-55* K from 


Kfre/D 


from O.67 to 16.67. In figure 3> the negative 


rk rk 

values of —2 2 as a function of are plotted for a similar range 

Kf r e /D 

of variables. This method of plotting has been used to avoid the diffi- 
culty of attempting to plot an Mg = 1 curve on a figure having Mq as 

an abscissa. In such a plot, the Mg = 1 curve would be a function of 
Kj>r e /D. As plotted in figure 3, the M e = 1 curve is independent of 
Kfre/D. 


In order to use figure 3 with only Mq known, it is necessary to 
find, by trial, the value of Mg that agrees with the known value of Mq. 

2 


A linear interpolation between given values of 
accuracy in the use of the charts . 




Kfrg/D 


is of sufficient 


Effect of K^r e /D 


As mentioned previously, the curves of figures 2 and 3 have been 
computed for Kpr e /D = 0.06. As a check on the effect of Kpr e /D, 
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calculations have been made for values of Kpr e /D of 0.02, 0.10, 0.24, 

and 0-3* These check points are shown in figures 2 and 3» At 

Kfr e jli = 0.10, the difference in exit Mach number between a calculated 


curve and the check point at the same 


M. t £ 


never exceeded 0.01 for 


Kf*e/ D 

any of the check points calculated. This discrepancy is well within the 
accuracy of the calculations. It appears, therefore, as though the 
charts of figures 2 and 3 are substantially independent of KfT e jI). As 

T 

further evidence of this independence, calculations at = 0.63, 

Mrp 2 T t,o 

— = 2.9, and Kfr e /D = 0.5 give a value of 

Kfrg/D ' 


M 0 = 0.24, k = 0.3, 


0 o - 0 e 

; — = 2.62 which agrees with the value obtained from figure 2(d) by 

Kfr e /D 

interpolation. 


Effect of Duct-Diameter Variation 

The charts may be used to compute the Mach number variation of a 
flow passage composed of a numiber of ducts having different diameters. 

In this case, the value of r e is the distance from the center of 
rotation to the exit of the duct being considered and Mq is based on 
the diameter of this duct. An application of the charts to this type 
of duct is given In the sample calculation. 

The preceding method can be applied to a tapering duct. However, 
if the duct tapers rapidly, the increments of duct length that must be 
used may approach the increments that would be used in Integrating 
equation (8b). In this case, there may be little or no saving in calcu- 
lation time between the two methods. 


Sample Calculation 

The application of the charts may best be illustrated by a sample 
calculation. Assume a duct consisting of two lengths of constant- 
diameter ducting. The first length extends from the 15-percent to the 
70-percent spanwiBe station and has a value of Kfr e /D = 0.10. The 
second length extends from the 70-percent to the 100-percent spanwise 
station and has a value of Kj=r e /D = 0.12. Assume that the ratio of 

the duct area of the second length to that of the first is 0.8l. The 
duct-tip Mach number Is assumed to be 0-9 and the duct flow is assumed 
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to be adiabatic (k = 0). Finally, assume that, at x = 0.15, the duct 
Mach number is 0.40 and the ratio of ambient temperature to toted tem- 
perature is O.65. 


Then, from equation ( 4b ) , 
x = 0.15, 


^t,o 

T t,i 


T 

O.997 or — £_ « O.65. At 
T t,o 


= [(0.15X0.9)] 2 (Q-7) 

%r e /DL i 5 (0.10)(0.15) 


0.850 


In order to calculate Mq for the first length of duct, assume that 
Mg = M^ = 0.40. Then, from figures 2(d) and 3(h), 


i*o - i*e 
K f r e /D 


-0.25 


and, by use of figure 1, 


Mq = 0.395 


For the first length of duct, the Mach number at x = 0.70 may be 
determined from figure 2(d) at 


Then, 


and 


' % 2 | 
Kf'Vv 0.70 


_(0-7)(0.9 )] 2 

0.10 


5.97 


0o - 0e 

W D 


2.45 


Mg = 0.345 

In the second length at x = 0.70, the Mach number Is 0.455 if a 
2-percent total -pressure loss is assumed at the contraction. In order 
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to calculate Mq for the second duct, assume that the Mach n umb er at the 
inlet to the second duct is an exit Mach number] that is. Mg = = 0 . ^55 • 

Thus , 

( \ = C(0-T)(0.9)] 2 = u 7 , 

\Kf r e/7 0 . 70 (0.12)(0.7) 

and 


T. 


t,o 


0.65 


Using figure 2(d) and assuming values of Mq until the correct value 
of M e is obtained yields Mq = 0.51* 

Then, at x = 1.0, 


From figure 2(d), 


and 


V 


JC^e/Dy 


0.81 

0.12 


= 6.75 


0o ~ 0e 
Kf r e/ D 


= 3.87 


M e = 0.38 

Since the mass flow and total temperature at the duct exit are known, 
any other flow property may he obtained. 


CONCLUDING REMARKS 


Charts have been developed for the analysis of flow in a whirling 
constant -diameter duct. These charts were presented for a range of 
variables which include the duct-inlet Mach number, duct-tip Mach number, 
friction factor, ratio of ambient temperature to duct-flow total temper- 
ature, ratio of duct-tip radius to duct hydraulic diameter, and a 
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dimensionless heat-transfer term. Hie method of using the charts is 
illustrated through the computation of a sample problem. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 3; 1957* 
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Figure 2.- Continued. 
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Figure 2.- Continued, 
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Figure 2.- Concluded. 





